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Fermi GBM as a Transient 
Monitor 


Colleen A. Wilson-Hodge 
Fermi GBM PI 
NASA/MSFC 


Outline 


GBM Instrument 

Brief transients (milliseconds to minutes) 

— On-board triggers and ground searches 

— Science of EM counterparts to GWs 

— GBM X-ray burst detections 

Long transients & variable sources (hours to 
years) 

— Earth occultation and accreting pulsar monitoring 
— Observational results 

Conclusions/Future work 


Fermi Gamma-ray Burst Monitor (GBM) 


r = Nal 


°GBM detectors are scintillating crystals 
with attached photomultipliers 


¢ Nal: 8 — 1000 keV 
¢BGO: 200 keV — 40 MeV 


Meegan et al 2009, Apj, 702, 791 


Fermi Gamma-ray Burst Monitor (GBM) 


a GBM Field of View 


Par 
ns { Fermi GBM provides gamma-ray 
SUlie\.e |context observations in the Multi- 
en 5 


Air 


= 87% uptime (off due to SAA), 

=» Views 69% of sky (Earth blockage), 

=" Views a particular point on the sky 
60% of the time, on average. 


GBM Localization Method 


High Expected Rates 


Low Rates 
Spacecraft Blockage 


¢ Localization is performed by comparing the relative observed rates 
from the GRB in each detector to the expected rates on a 1 degree 
grid 

¢ This requires an assumption of the spectrum, and the sky grid limits 
to a statistical minimum uncertainty of 1 degree radius 


GBM Detector / Instrument Response 


GRB Photon Spectrum Instrument Response Instrument Background Observed Data 


GBM Data 


Data Type Energy Resolution 
TRIGDAT 1024/256/64 ms 


CSPEC 4096/1024 ms 128 channels 
cTTE (New 


TRIGDAT: used primarily for localization & quick look (triggers only) 


CTIME: temporal analysis 

CPSEC: spectral analysis 

Initially TTE was available ~30s pre-trigger - ~300 s post-trigger 
Continuous TTE (CTTE) implemented on November 26, 2012 


Brief Transients 
(milliseconds to minutes) 


GBM On-Board Triggering 


¢ GBM triggers when 2 or more detectors exceed background by n 
sigma over timescale t in energy band e. 


¢ 70 algorithms operating simultaneously. 
—- 45<n</7.5 
— 16ms<t<8.096s 
— e=oneof 25-50 keV, 50 - 300 keV, 100 - 300 keV, > 300 keV 


= What does GBM trigger on? 


A. von Kienlin MPE 
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748 TGFs 


Quarter 


446 Other, 883 particles 


1126 Solar Flares 


270 SGR bursts 


GBM Triggered GRBs 


2000 Fermi GBM GRBs 


13% Also seen by Swift 
6% Also seen with LAT 


2000 GBM GRBs 
266 Swift GRBs 


GBM GRBs 
121 LAT GRBs oe 


52% in LAT FoV 
48% Outside LAT FoV 


2000 GBM GRBs 
121 LAT GRBs 
In Field-of-view of LAT (1045) 
Out of Field-of-view of LAT (955) 


Short and Long GRBs 


2000 Fermi GBM GRBs 
Wes - Drege: 9 - -- 


Each year GBM sees: weer 
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Short GRB / CBC Association 


Jet-ISM Shock (Afterglow) 
Opiical (hoars—days) 


Radio (weeks—years} 


Ejecta-ISM Shock 
Radio (wears) 


Short GRB / CBC Association 


Short gamma-ray burst 
Coe iclaelalecmelll cities) 
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Stars* in GW 


a compact 


binary system fea * In-spiral confirms CBC progenitor model 

2egIN CO Spita 1 . 

inward. ¢ Information about binary system parameters 
* precise merger time 


¢ standard candle -> luminosity distance 


eventually 
rae) (ellale p 


EM 
Vis ouking wis ¢ EM energetics 
Lorin yaa ¢ X-ray or optical afterglow gives precise location 
sa hs Breaks degeneracy in binary parameter estimation 
Host galaxy/redshift 
Local environment information 
With many: jet opening angle 


*Possibly neutron stars. 


Joint Sub-threshold Searches 


Ideal Scenario Bright GBM Bright LIGO 


GW150914 Sub-threshold Bright LIGO 
Scenario GBM 


Typical more —_ Bright GBM Sub-threshold 
distant short LIGO 

GRB 

Both Sources Sub-threshold Sub-threshold 
Faint GBM LIGO 


The GBM and LIGO teams have been working together to develop automated pipelines 
to search for sub-threshold signals 


In all cases, the presence of a signal in GBM or LIGO, can raise the significance of the 
signal being real in the other instrument. 


A confident gamma-ray signal enables detecting a fainter gravitational wave signal, 
increasing the GW detection distance limit, in turn increasing the event rate by a factor 
of the distance cubed 


GBM Ground Searches 


¢ Un-targeted search 

— Looks for signals too faint to trigger on-board 

— Needs no input from other instruments 
¢ Targeted search 

— Seeded with a time and optional sky map 

— Inputs generally from LIGO, neutrino detections, etc. 
¢ X-ray Burst search 


— Manual search, resulting from data cleaning for pulsar 
monitoring 


Un-targeted Search (1) 


httos://gammaray.nsstc.nasa.gov/gbm/science/serb se 
arch.html 


Developed to search for faint short GRBs 


Uses CTTE data, 2us time resolution, 128 energy 
channels 


18 timescales: 64 ms to 32s 
— On-board: 16 ms to 8.096 s 
Time series are made 4 times, offset in phase 
5 energy ranges (optimized on triggered sGRBs) 


Fit a background using cubic splines and filter out bad 
background fits 


Fast, efficient, runs on hourly CTTE data as it arrives 


Un-targeted Search (2) 


= Test for statistically significant excesses in two 
Nal detectors: 
o Require 2 2.50 in the best Nal detector & 
2 1.250 in the second best Nal detector, 
o Require Poisson probability < 1E-6, including 
trials factor for Nbins in one day, 
o Other trials factors not included. 
= Require the detector pair be valid for a distant 
point source. 


lg 


Un-Targeted Search Results 


Galactic coordinates Spacecraft Coordinates 


318 short, hard candidates (known triggers omitted) in 46 months. 
~80 per year, twice the rate of GBM triggered sGRBs. 


Verification 1/2 
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Positive candidates 
(blue) versus 
Negative deviations 
(red) 
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Un-Targeted Search Verification 


» ~1/4 have signals in more than 2 GBM detectors. 
= Swift GRB 140606A: not a GBM trigger, easily detected at P = 1E-20. 
= INTEGRAL Anti-Coincidence Shield (ACS): 

= ~1/3 of spectrally-hard candidates are detected by the ACS 

= GBM triggered sGRBs: ~50% detected with the ACS. 
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A very average candidate: signal in only 
two detectors & P=2E-7. 
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ACS - sGRB20170125a 
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Not detected in INTEGRAL ACS 


Targeted Search of GBM data to GW events 


Developed during LIGO S6 observing run (Blackburn et al. 2015, Goldstein et al. 2016) 


¢ Coherent search over all 14 GBM detectors (Nal and BGO) 


> 


> 
> 


seeded with time & (optionally) sky location of 


any LIGO/Virgo candidate event Coherent search over GBM detectors 
over user-specified time window (+30 s) imme 
fey 


estimate of background rate by polynomial to 
local data outside the foreground interval 


time 
GBM search window (~!m) 


@ For each template spectrum (soft, medium & hard) and sky location 


> 
> 
> 
> 


> 


Each model spectrum is folded through the detector response to determine detector counts 
Detector counts for each energy channel are weighted according to the modeled rate 
and inverse noise variance due to background 


Weighted counts from all Nal and BGO detectors are summed to obtain a signal-to-noise optimized 
light curve for that model 


Each model is assigned a likelihood by the targeted search based on the foreground counts 


@ Candidates are ranked by a Bayesian likelihood statistic 


@ Will reveal short-duration candidates between 0.256 s to 8.192 s (CTIME) 


GW150914-GBM 


Targeted search around GW150914: 
Initial 60s (+ 30s) search window (selected a priori) 


2 candidates 
° Soft transient: Ty + 115, 2s long: Gal.Cent. region 
¢ Hard transient: Toy +0.4s, 1s long: GW150914-GBM 


TITLE: GCN CIRCULAR 

NUMBER: 18339 

SUBJECT: LIGO/Virgo G184098: Fermi-GBM ground-based follow-up 
DATE: 15/09/20 01:46:08 GMT 

FROM: Lindy Blackburn at CfA <lindy.blackburn@ligo.org> 


Lindy Blackburn (CfA), Michael S. Briggs (UAH), Eric Burns (UAH), Jordan Camp 
(NASA/GSFC), Nelson Christensen (Carleton College), Valerie Connaughton 
(USRA), Adam Goldstein (NASA/MSFC), Tyson Littenberg (UAH), John Veitch 
(Birmingham), Judith Racusin (NASA/GSFC), Peter Shawhan (UMD), Leo Singer 
(NASA/GSFC), Binbin Zhang (UAH) 


We report on a sub-threshold targeted followup of LIGO candidate event G184098 
in Fermi-GBM survey data for bursts between 0.256s and 8s in duration, and 
covering a range of GRB spectral models. Although there was no on-board GBM 
trigger at the time of the event, Fermi-GBM was exposed to a large fraction of 
the LIGO sky position and thus we searched offline data for untriggered 
events. The GBM FOV is blocked by the Earth which occults 67 degrees from (RA, 
DEC) = (355.14, -21.23). Thus GBM observation is able to cover about 87.8% of 
the cWB sky posterior, and 91.5% of the LIB posterior. We scanned several 
minutes of GBM live-time centered on the GW event time using a pipeline 
developed specifically for following-up LIGO-Virgo events in GBM archival data 
during the LIGO-Virgo S6/VSR3 run [1]. 


The search identified a possible transient beginning at 150914 09:50:45.8, 
about 0.4s after the reported LIGO burst trigger time of 09:50:45.39, and it 
lasted for about 1 second. The intrinsic time resolution for this search was 
0.256s. Of the three GRB model spectra tested in the search, the event was 
best matched to the one corresponding to the hardest spectrum. Using GBM 


Raw count rates: 
Sum of all GBM detectors: 12 x Nal + 2 x BGO 
Nal: 50-980 keV / BGO: 420 keV — 4.7 MeV 
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Counts per Second 
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Model-dependent count rates: 
Raw count rates weighted & summed to max signal-to- 
noise for a modeled source 
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GW150914-GBM 


Targeted search around GW150914: 
— Initial 60s (+ 30s) search window (selected a priori) 


— 2 candidates 
° Soft transient: Tay + 115, 2s long: Gal.Cent. region 
* Hard transient: Toy + 0.45, 1s long: GW150914-GBM 
0.2% probability of occurring by chance (2.90) 


False Alarm Probability Calculation: 


: 1 = 27 hard events in 218821.1s of GBM live 
ane factor of 3 for spectra searched, 90% confidence 


-- between GW TO and 
GBM event start 


Effective ~ for bins/durations 
Searched: 
30 s: max offset (search window) 


~ to account for offset in 0.256 s: min CTIME bin 


time in either direction 


Raw count rates: 
Sum of all GBM detectors: 12 x Nal + 2 x BGO 
Nal: 50-980 keV / BGO: 420 keV — 4.7 MeV 


Counts per Second 
3 


6000 


0 
Seconds from GW TO 


Model-dependent count rates: 
Raw count rates weighted & summed to max signal-to- 
noise for a modeled source 


flux [counts/s] 


relative time [s] 


Characteristics of GW150914-GBM 


e Unusual detector pattern: 
nearly equal count rates in all Nal detectors 


— Localization: source direction underneath the spacecraft, 
163° to the spacecraft pointing direction 


Counts per Second 


800 
-10 -5 Oo 5 10 -10 -5 0 5 10 
Seconds from GW TO 


Nals: 


50 — 980 keV At nti 
BGOs: wa i 


420 keV — 4.7 MeV 


ne 


° 
Seconds from GW TO 


o deviation from a background fit 


Nal 0 Nall Nal2 Nal3 Nal4 Nal 5 
1.31 1.81 0.64 1.05 2.42 1.68 
Nal 6 Nal7 Nal8 Nal9 Nall0 Nal lil 
1.31 1.64 1.45 2.20 1.61 0.66 
BGOO BGO1 
2.25 2.56 


Characteristics of GW150914-GBM 


e Unusual detector pattern: 
nearly equal count rates in all Nal detectors 


— Localization: source direction underneath the spacecraft, 163° to 
the spacecraft pointing direction 

— If association with GW150914 was true: shrink LIGO localization 
by 2/3 
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Characteristics of GW150914-GBM 


e Unusual detector pattern: 


-1.0 


nearly equal count rates in all Nal detectors 


¢ Energy spectrum: 


Localization: source direction underneath the spacecraft, “1.2 
163° to the spacecraft pointing direction 

If association with GW150914 was true: shrink LIGO 
localization by 2/3 


Power Law (Photon) Index 


Peaking in BGO energy range 


Best fit simple PL with index -1.4 (average for sGRBs), Ee 
Fluence 2.4 x 107 erg cm? (weaker than average for 
sGRBs) 
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@ All Short GRBs 
®@ Best fit with PL 
@ Gw150914-Gem 
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Fluence [10-1000 keV] (erg cm’) 


Times: 0.384: 1.408 s 
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Association with GW150914? 


¢ Evidence for e Evidence against: 

— 3sigma False Alarm Probability — Low significance 

— GBM signal localized to a region — Lack of corroboration by other 
consistent with the LIGO sky map experiments 

— Cannot be attributed to other — Nature of the LIGO event is a BH- 
known astrophysical, solar, BH merger 
terrestrial or ial ae? 
activity 


Fermi Orbit 
arise aes ees 


Lightning 
(TGFs/TEBs) 


Galactic 
Sources 


Magneto- 5 9 
spheric . : 


The most likely explanation is a short GRB ... 
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GBM Observations of GW Events 


GW 151226 
(Abbot et al. 2016b) 

¢ BH+BH Merger 
° 144&75Me 

¢ 440 Mpc 


GW 150914 LVT 151012 
(Abbot et al. 2016a) (Abbot et al. 2016a) 

¢ BH+BH Merger * Candidate BH+BH 
* 36&29Me * 23&13M6 

¢ 410 Mpc ¢ 1100 Mpc 


= 
— 
Counts per Second 


. -30 -20 -10 0 10 #420 #30 
= Seconds from tyyr 


0 5 
Seconds from GW TO 


=30 -20 -10 0 10 #20 #30 
Seconds from tew 


(Connaughton et al 2016) (Racusin et al 2016) 


GW150914-GBM, a 2.90 event consistent with a short GRB 


— Not predicted by theoretical models 


No gamma-ray detections for LVT151012 or GW151226 — not constraining 


— 32% and 17% of LIGO localization region blocked by Earth for GBM 
— Backgrounds were 18% and 3% higher in GBM 
— Distance for L\VT151012 was 3x larger 


— lf gamma-ray emission is in a jet, only 15-30% would be pointed toward Earth 


Need more events before we can say more! 


Type 1 X-ray Burst 


Neutron star accreting matter from a low mass companion at low mass accretion rate. 


¢ Three types 


— Normal 
¢ 10-100 s 
° H, He 
¢ L~ 10°9-— 10% erg/s 
— Long 
¢ 10-30 Minutes 
owe 
¢ L~ 10* erg/s 
— Super 
¢ Hours — Days 


: gE 
Mw 2% Eddington ¢ L~ 10% erg/s 


GBM X-ray Burst Search 


Visual Inspection of CTIME Data 


12 Nal detectors GBM Pulsar Project 


12-25 keV PI C. Wilson-Hodge 
8 second bins 


CTIME Editor 


4U 0614+09 


Initiated March 12, 2010 


Identification Process 
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Spectral 


Location 


What does GBM see? 


4U 0614+09 es ener 
10° 7 
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10° _Photospheric Radius 
Expansion (PRE) Touchdown 


Cape: %. 


F (10° erg/cm?/s) i 


kT (keV) 


ASM PCA 


Kuulkers et al. 2010 © 


R(km) 
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Three Year X-ray Burst PAAUE: 


——+—__+—__+--y 


https://gammaray.nsstc.nasa.gov/gbm/science 


/xrb. html 


Jenke et al. 2016, ApJ, 826, 228 


752 Thermonuclear XRBs 
267 Transient Events from accretion flares 
65 Untriggered GRBs 


GBM is sensitive to photospheric radius 
expansion (PRE) bursts 


1.4 PRE bursts per day within 10 kpc 
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Associations for Low M,., Accretors 


Associations For Low M,,, Accretors 


@ Known Type | X-ray Bursters 


33 - 4U 0614+09 eee vig ok Targeted Type 1 X-ray Bursters 


4 - SAX J1818.74+1424 
2 - UW Crb 

2 —IGR J17062-6143 
1 — XB 1940-04 

1 —Ser X-1 

1 — MAXI J1421-613 
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Locations are poor. Must use 
MAXI rates to determine if 
potential source is active. 

Automatic checking. 


WN 


Long Transients 
(hours to years} 


Accretion Powered Pulsar Monitoring 


¢ Blind search 

— For unknown sources and unmonitored transients 
¢ Dedicated search 

— Search around known frequencies 

— Currently monitoring 39 systems (36 detected) 
¢ GBM Advantage 

— Typically > 40,000 s of on-source time per day! 
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Roche lobe overflow : 
Wind accretion Accretion from a Be star’s 


circumstellar disk 


httos://gammaray.nsstc.nasa.gov/gbm/science/pulsars.html 


GRO J1744-28 
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MAXI J1409-619 


55000 55500 56000 56500 57000 57500 
Time (MJD) 


EXO 2030+375 (P.=42 s,P,,=46 d) 


Discovered during a giant outburst in 1985 with 
EXOSAT; Second giant outburst in 2006 (RXTE) 


Transitioned to spin down in 1995 and again in 
2016 

Abruptly shifted in outburst orbital phase in 1995 
accompanied by a drop in outburst flux; again in 
2016 

Detected outburst at nearly every periastron 
passage since 1991, unlike most Be X-ray binaries 
Correlated peak flux and orbital phase of 
outburst peak — delay of accretion from Be disk 
onto NS accretion disk? 


orb 


EXO 2030+375 Torque Reversal 
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Spin Frequency (mHz 
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12-50 keV Pulsed Flux 
(keV cm* s"') 
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oe in Be X-ray binary EXO 2030+375 


EXO 2030+375 April 2016 
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GBM Discovers rare torque reversal 
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EXO 2030+375 Long Term Behavior 
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Orbit Determination — Flux based torque model 


Torque Model |; 


Frequency (cycles/day) 


v = frequency derivative 


B2/7 79/7 Bip = Magnetic Field 
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Orbit Determination 
Polynomial Torque model 
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‘Occultation steps’ occur in the detector count rates as sources 
rise above and set below the Earth’s horizon... 


)) Crab Nebula — 
detected to 300 keV 


Occultation 
step 


C. A. Wilson-Hodge, P. Jenke, et al., “Three years of Fermi GBM Earth 
Occultation Monitoring: Observations of Hard X-ray/Soft Gamma-Ray 
S1o]0] (6(-\- ela OU bo -40 GIO (408 2) 


GBM Observations of V404 Cyg 


10 Meoblack hole only 2.4 kpc away 

Last X-ray outburst observed with Ginga in 1989 
Two confirmed optical outbursts in 1938 and 
1956 

169 GBM triggers over 13 days starting June 15, 
2013 

73 distinct flaring episodes 

Reached 30 crab with emission up to 300 keV 


GBM observed the entire outburst. It can be 
analyzed using both triggered data and Earth 
occultation 
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V404 Cyg is predominately 
in the Low/Hard State 
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Flux 12-300 keV (Crab) 


RXTE PCA rate (c/s/PCU) 


Haninoss 


HR = 8.6-18 keV/5.0-8.6 keV 


HR: 12-—25/8-12 keV 


Long-term Hard X-ray Variations in the 
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Searches for multi-messenger 
counterparts using Earth occultation 


Using the Earth occultation technique, we search for 
new sources by measuring fluxes for source positions 
covering the LIGO arc or for the neutrino detection 
postion. This is a search for extended (time) emission. 
We performed these searches for GW150914 
(Connaughton et al 2016), GW151226, and LVT151012 
(Racusin et al. 2016) using + 1 day, 1 month, and 1 year 
if data were available. 

We have performed similar searches for neutrinos 


To date we have not detected any new sources in these 
searches 


Summary 


GBM detects a broad range of transients, on 
timescales from milliseconds (TGFs) to years (X- 
ray binaries and Crab Nebula variations) 

GBM is especially well suited to detecting short 
GRBs due to its wide field of view, high duty cycle, 
and broad energy range. 

GBM is the most prolific detector of short GRBs 
available and we are eagerly awaiting detections 
of gravitational waves from NS-NS mergers! 

In the meantime, there is plenty of science to do 
with GBM data. 


